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: [8 i t RNA ( circRNA ) 在 可 变 剪 接 、 转 录 调 控 和 来 源 基因 的 表达 调控 等 方面 具有 重要 
能 。 本 研究 分 析 意 大 利 蜜蜂 Apis mellifera ligustica 工蜂 中 肠 circRNA 的 数量 种类、 结构 特 
征 和 作用 ， po 分 析 circRNA 的 调控 网 络 探索 circRNA 的 调控 功能 。【 方 法 】 在 实验 室 条 
件 下 人 工 饲 养 意大利 蜜蜂 工蜂 ,利用 circRNA-seq 技术 对 意大利 蜜蜂 7 Fo 10 日 龄 成 年 工蜂 中 肠 样 
品 进行 深度 测序 。 利 用 find_cire 软件 从 质 控 后 的 数据 中 预测 circRNA, iit BLAST R3] GO 和 
KEGG 数据 库 , 对 circRNA 的 来 源 基 因 进 行 功能 和 代谢 通路 注释 。 利 用 TargetFinder. 软件 预测 
circRNA 靶 向 结合 的 miRNA 及 miRNA 靶 向 结合 的 mRNA ,通过 Cytoscape v.3.2. 1 软件 对 circRNA- 
miRNA 和 circRNA-miRNA-mRNA 调控 网 络 进行 构建 及 可 视 化 。 通 过 设计 背靠背 引物 和 线性 扩 增 
引物 RT-PCR 对 预测 出 的 circRNA 进行 验证 。【 结果 ] 意 大 利 蜜蜂 工蜂 中 肠 样品 的 测序 平均 得 到 
136 463 071 条 clean reads, 去 除 rRNA 后 各 样品 的 anchor reads 均 在 136 779 122 条 及 以 上 。 共 预测 
出 10 833 个 circRNA ,长 度 主要 介 于 15 ~1 000 nt; E38 circRNA 的 类 型 丰富 ,其 中 已 注释 的 外 显 子 
circRNA 数量 最 多 , 分布 在 西方 蜜蜂 1 号 染色 体 的 circRNA 数量 最 多 ,其 次 为 8 号 染色 体 。 
CircRNA 的 来 源 基 因 可 注释 到 包括 结合 、 细 胞 进程 和 细胞 在 内 的 45 个 GO 条 目 ,以 及 包括 内 吞 作 
、 内 质 网 蛋白 加 工 及 核糖 体 在 内 的 121 条 KEGG 代谢 通路 ,表明 circRNA 在 意大利 蜜蜂 工蜂 中 肠 
i EK 发育 新陈代谢 和 细胞 生命 活动 等 生物 学 过 程 中 发 挥 重 要 作用 。 进 一 步 构 建 circRNA- 
miRNA 和 circRNA-miRNA-mRNA 调控 网 络 , 分 析 结 果 显 示 部 分 circRNA 可 能 作为 竞争 性 内 源 RNA 
吸附 结合 microRNA ,从 而 调控 基因 的 表达 水 平 。 最 后 ,对 随机 选择 的 3 个 circRNA 的 RT-PCR 结果 
aa ian 【结论 】 本 研究 对 意大利 蜜蜂 工蜂 中 肠 中 的 circRNA 进行 预测 、 分 析 及 鉴定 。 
究 结果 提供 了 中 肠 circRNA 的 数量 、 种 类 、 结 构 特 征 、 作 用 和 调控 网 络 的 信息 ,揭示 了 circRNA 可 
i. 于 来 源 基 因 和 作为 竞争 性 内 源 RNA 在 意大利 蜜蜂 工蜂 中 肠 的 生长 发 育 和 免疫 防御 中 
发 挥 作用 ,为 深入 研究 circRNA ek 蜜蜂 中 肠 发 育 及 胁迫 响应 过 程 中 的 功能 黄 定 了 基础 。 
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Abstract: [ ^im] Circular RNA (circRNA) plays a key role in alternative splicing, transcription 
regulation and expression regulation of source genes. This study aims to analyze the quantity, variety, 
structural characteristics and function of circRNAs in the midgut of Apis mellifera ligustica workers, and to 
explore the regulatory function of circRNAs via constructing and analyzing regulatory networks. 
[Methods] A. m. ligustica workers were reared under laboratory conditions, and the midgut samples 
from 7- and 10-day-old workers were subjected to deep sequencing using circRNA-seq technology. 
CircRNAs were predicted from sequencing data after quality control using find, cire software. Source 
genes of these circRNAs were annotated to GO and KEGG databases to gain function and pathway 
annotations via BLAST. The target mRNAs of circRNAs and miRNAs were predicted with TargetFinder 
software, and the regulatory networks between circRNAs and miRNAs and between circRNAs, miRNAs 
and mRNAs were constructed and visualized using Cytoscape v. 3. 2. 1 software. The predicted circRNAs 
were validated by RT-PCR with the divergent and convergent primers designed. [Results] Sequencing of 
midgut samples from A. m. ligustica workers produced a mean of 136 463 071 clean reads, and for each 
sample over 136 779 122 anchor reads were obtained after removing rRNA. A total of 10 833 circRNAs 
were predicted, and their length ranged mainly from 15 to 1 000 nt. The types of these circRNAs were 
abundant, and the largest one was annotated to be exonic circRNA. The number of circRNAs distributed 
on chromosome 1 of Apis mellifera was the most and that distributed on chromosome 8 was the second. 
The source genes of circRNAs could be annotated to 45 GO terms including binding, cellular process, 
and cell, as well as 121 KEGG metabolic pathways including endocytosis, protein processing in 
endoplasmic reticulum, and ribosome, suggesting that circRNAs may play key roles in such biological 
processes as growth, development, metabolism, and cellular activity in the midgut of A. m. ligustica 
workers. Furthermore, the regulatory networks between circRNA and miRNA and between circRNA, 
miRNA and mRNA were constructed, and the analysis result demonstrated that partial circRNAs could 
bind microRNAs as competitive endogenous RNAs (ceRNAs). Finally, the true existence of randomly 
selected three circRNAs was validated by RT-PCR. [Conclusion] In this study, we predicted, analyzed, 
and verified circRNAs in the midgut of A. m. ligustica workers. Our findings provide the data of the 
quantity, variety, structural characteristics, roles, and regulatory network of circRNAs in the midgut, 
indicating that circRNAs can play roles in the developmental and immune defense processes in the midgut 
of A. m. ligustica workers through affecting source genes and acting as ceRNAs. This study lays a 
foundation for further studying the roles of circRNAs in the development and stress response of the midgut 
of A. m. ligustica. 
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意大利 蜜蜂 Apis mellifera ligustica ( 简称 意 蜂 ) 
属于 西方 蜜蜂 Apis mellifera ,具有 群 势 强 .产量 高 不 
喜 分 蜂 等 优良 生产 性 能 ,因而 被 广泛 用 于 世界 各 国 
的 养 蜂 生产 (Jensen et al., 2010) 。 早 在 2006 年 , 西 
方 蜜蜂 的 基因 组 信息 就 已 完成 测序 并 公布 
(Honeybee Genome Sequencing Consortium, 2006 ) , 
此 后 西方 蜜蜂 的 分 子 生物 学 和 组 学 研究 获得 了 长 足 
发 展 。 

非 编 码 RNA (non-coding RNA, ncRNA) 起初 被 





























认为 是 基因 组 的 “暗物质 ”和 基因 转录 的 “噪音 ”。 
近年 来 , 随 着 测序 技术 的 飞快 发 展 和 相关 研究 的 不 
WRA ncRNA 已 成 为 备 受 关注 的 研究 热点 。 根 据 
长 度 和 形状 ,ncRNA 可 进一步 分 为 小 RNA .长 链 非 
编码 RNA (long non-coding RNA, IncRNA) 和 环 状 
RNA(circular RNA, circRNA) f, CircRNA 是 新 近 
发 现 的 一 类 ncRNA ,通过 外 显 子 或 (和 ) 内 含 子 的 反 
向 剪 切 形成 环 化 RNA 分 子 ( Meng et al., 2017) ,在 
不 同 物种 中 具有 保守 性 .丰富 性 稳定 性 和 组 织 表 达 
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特异 性 等 特点 (He et al., 2017) ,已 被 证 明 具 有 多 种 
生物 学 调控 功能 ,包括 作为 微小 RNA ( microRNA, 
miRNA ) 海 绵 ,形成 RNA- 蛋 白质 复合 物 ,以 及 调控 
靶 基 因 的 转录 和 可 变 剪 接 等 (Qu et al., 2015) 。 较 
之 线性 RNA,circRNA 缺少 5" 帽子 和 3"' 尾 巴结 构 , 
而 稳定 性 更 高 且 能 抵抗 RNase R 的 消化 ,成 为 理想 
的 生物 标志 物 (Perkel, 2013)。 借 助 高 通 量 测序 技 
术 和 生物 信息 学 分 析 方 法 ,迄今 已 在 人 类 Homo 
sapiens ( Jeck et al., 2013; Salzman et al., 2013; 
Memezak et al., 2013) , 3. 44 Danio rerio ( Shen et 
al., 2017) 水稻 Oryza sativa (Lu et al., 2015) .蜜蜂 
PRHE fal Ascosphaera apis (Guo et al., 2018) 秀丽 隐 
fF 2 E Caenorhabditis elegans ( Liu and Chen, 
2015) . 古 细 菌 (Danan et al., 2012) 等 动 植物 和 微 生 
物 中 预测 和 鉴定 出 大 量 的 circRNA ,但 昆虫 的 
circRNA 研究 目前 尚 处 于 初级 阶段 。Westholm 等 
(2014 ) 对 Hz ow Drosophila melanogaster 的 
circRNA 进行 全 基因 组 分 析 , 发 现 果 蝇 的 circRNA 
具有 组 织 表达 特异 性 ,并 且 可 作为 果 蝇 衰老 的 潜在 
生物 标志 物 。Gan 等 (2017 ) 3 et Xf KA Bombyx 
mori 中 肠 和 后 肠 进行 高 通 量 测序 预测 出 来 源 于 
1 727 个 基因 的 3 155 个 circRNA ,进一步 分 析 发 现 中 
肠 和 后 肠 中 的 来 源 基 因 的 功能 和 代谢 通路 
(pathway) 注释 信息 类 似 。Chen 等 (2018 ) 对 意 蜂 蜂 
王 的 卵 梨 组 织 中 的 cireRNA 进行 预测 和 分 析 , 发 掘 
出 12 211 个 circRNA ,并 通过 对 刚 产 卵 蜂王 、 处 女 
王 限制 产 卵 蜂王 和 产 卵 恢复 蜂王 进行 比较 分 析 推 
测 差异 表达 的 circRNA ( DEcircRNA ) 可 通过 竞争 性 
结合 miRNA 在 蜂王 卵巢 组 织 活化 及 产 孵 过程 中 发 
挥 作用 。 

笔者 所 在 课题 组 前 期 在 mRNA 组 学 水 平 对 意 
蜂 幼 虫 肠 道 的 发 育 机 理 、 胁 迫 应 答 机 制 进 行 了 一 系 
列 研究 ,提供 了 意 蜂 幼虫 肠 道 发 育 过 程 中 的 基因 表 
达 谱 和 差异 表达 信息 ( 郭 寄 等, 2018b ) ,解析 了 意 蜂 
幼虫 肠 道 啊 应 球 吉 于 胁 迫 的 免疫 应 答 ( Chen et al., 
2017) ;此 外 还 在 IncRNA 组 学 水 平 解析 了 意 蜂 工蜂 
中 上 肠 发 育 过 程 中 的 IncRNA 25 5e xe ^ i Cp EARS, 
2018a)。 日 前 ,蜜蜂 circRNA 的 相关 研究 十 分 荐 乏 ， 
AREE circRNA 的 结构 和 功能 的 信息 极为 有 限 。 
本 研究 结合 circRNA-seq 技术 和 生物 信息 学 分 析 方 
法 对 意 蜂 工蜂 中 肠 的 circRNA 进行 系统 的 预测 和 结 
构 分 析 , 通过 来 源 基 因 分 析 、 调 控 网 络 分 析 进 行 
circRNA 的 功能 预测 ,并 通过 设计 特异 性 的 背靠背 
引物 (divergent primer) 和 线性 扩 增 引物 (convergent 
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primer) 对 预测 出 的 cireRNA 进行 验证 。 研 究 结果 
可 提供 意 蜂 工蜂 中 肠 circRNA 的 数量 .种 类 结构 和 
作用 的 重要 信息 ,为 蜜蜂 的 ncRNA 信息 提供 有 益 补 
充 , 同 时 也 为 深入 研究 circRNA 在 中 肠 发 育 及 胁迫 
响应 过 程 的 作用 建立 基础 。 




















1 材料 与 方法 


1.1 供 试 昆虫 

意大利 蜜蜂 工蜂 取 自 福建 农林 大 学 蜂 学 学 院 教 
学 蜂 场 。 从 外 观 健 康 且 群 势 较 强 的 蜂 群 中 选择 蜂 子 
状况 良好 的 成 熟 封 盖子 脾 , 迅 速 提 至 实验 室 , 放 入 
34( +0. 5) C 培养 箱 , 将 刚 出 房 的 工蜂 ( 记 为 0 日 
龄 ) 放 入 四 周 打 孔 的 干净 塑料 盒 (10 AE) ,每 个 逆 
料 盒 上 方 搬 入 一 支 逆 有 50% (w/v) Fo ba BK B5] S 
器 。34 +0. 5C 饲养 工蜂 至 10 日 龄 。 每 日 检查 工蜂 
存活 情况 ,及 时 清理 死亡 工蜂 。 本 实验 进行 3 次 生 
物 学 重复 。 待 工蜂 7 HET 10 日 龄 时 ,在 超 净 台 用 
干净 杀 子 拉 取 工蜂 中 肠 , 放 和 人 RNA-Free 的 EP 管 ， 
经 液 氮 速冻 后 迅速 转移 至 -80% 超低温 冰箱 保存 备 
JH. 7 日 龄 成 年 工蜂 中 肠 样品 的 3 个 重复 : Am7-1 ， 
Am7-2 和 Am7-3; 10 日 龄 工蜂 中 肠 样品 的 3 个 重 
fi: Am10-1, Am10-2 和 Am10-3 。 
1.2 RNA 提取 .cDNA 建 库 及 circRNA-seq 

(1) 利 用 RNA 抽 提 试剂 盒 AxyPrep™ Multisource 
Total RNA Miniprep Kit (TaKaRa 公司 ,日 本 ) 抽 提 
1.1 节 蜜 蜂 中 肠 样品 的 总 RNA ,为 最 大 限度 地 保留 
所 有 非 编 码 RNA(ncRNA ) ,去 除 核 糖 体 RNA 后 的 
mRNA 和 ncRNA 用 裂解 缓冲 液 随 机 打 断 为 小 片段 ， 
作为 模板 用 六 碱 基 随 机 引物 、 缓 冲 液 dNTPs, RNase 
H 和 DNA polymerase I 合成 cDNA 第 2 链 ;(2) 经 过 
QiaQuick PCR 试剂 盒 (Qiagen 公司 ,德国 ) 纯化 并 加 
EB 缓冲 液 洗 脱 经 末端 修复 .加 碱 基 A, 加 测序 接头 ， 
然后 通过 尿 喀 啶 -N- 糖 基 化 酶 (UNG ) 降解 cDNA 第 
2 链 ;(3) 消 化 产物 经 琼脂 糖 凝 肌 电 泳 和 PCR 扩 增 ， 
PCR 产物 经 高 碘 酸 钠 ( Sigma 公司 ,美国 ) 处 理 ; (4) 
委托 广州 基 迪 奥 生物 科技 有 限 公 司 对 上 述 cDNA X 
库 进行 双 端 ( paired-end ) 测序 ,测序 平台 为 Humina 
HiSeq 4000。 测 序数 据 已 上 传 NCBI SRA 数据 库 ， 
BioProject 号 : PRJNA406998 。 
1.3 ”测序 数据 的 质 控 、circRNA 的 预测 及 结构 特征 
分 析 

对 于 测序 得 到 的 raw reads, 去 除 未 知 碱 基 (N) 
的 比例 大 于 1096 的 、 质 量 值 0<10 的 和 碱 基数 占 整 
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条 读 段 数 的 5096 以 上 的 读 段 ,过 滤 得 到 的 clean 
reads 用 于 后 续 的 数据 分 析 。 

利用 TopHat 软件 (Kim et al., 2013) 将 各 中 肠 
样品 的 clean reads 与 西方 蜜蜂 参考 基因 组 (Amel_ 
4. 0) (Honeybee Genome Sequencing Consortium, 
2006) 进行 比 对 , 从 比 对 结果 中 提取 unmapped 
reads ,然后 截取 每 一 条 unmapped reads 的 两 端 ( 默 
认 20 bp) ,得 到 anchor reads。 继 而 用 短 序 列 比 对 软 
{F Bowtie2 (Langdon, 2015) 将 anchor reads 比 对 参 
考 基 因 组 ,将 得 到 的 比 对 结果 提交 给 find. cire 软件 
(Memezak et al., 2013) ,从 而 对 circRNA 进行 筛选 
和 鉴定 ,筛选 条 件 包 括 : (1) breakpoints = 1, BR 
留 有 且 只 有 1 个 清晰 breakpoint 的 环 状 RNA; (2) 
anchor overlap <2 , 即 每 条 reads 的 两 个 anchor reads 
比 对 到 基因 组 上 的 位 置 overlap 不 能 超过 2 bp; (3) 
edit 和 2 , 即 只 允许 2 bp 错 配 ;(4)n uniq >2, 即 uniq 
reads 大 于 2 4;(5) best qual A >35 BK best qual B > 
35 , 即 每 条 reads 的 其 中 一 个 anchor reads 比 对 到 基 
因 组 上 最 好 的 mapping 结果 要 比 其 排 第 二 的 结果 的 
分 值 高 35 分 以 上 ;(6)n unig > int( samples/2) , Bl xc 
持 该 circRNA 的 uniq reads 要 大 于 总 样品 数 的 二 分 
之 一 ;(7)circRNA 的 长 度 小 于 100 kb. 

对 circRNA 的 长 度 进行 统计 ,绘制 cireRNA 的 
长 度 分 布 图 。 根 据 find. circ 给 出 的 circRNA 信息 和 
参考 基因 组 的 转录 组 注释 ,对 鉴定 到 的 circRNA 进 
行 环 化 类 型 统计 ,包括 (1) 已 注释 的 外 显 子 , 即 
circRNA 的 breakpoint 位 于 已 知 基 因 的 一 个 外 显 子 
的 起 始 和 另 一 个 外 显 子 的 终止 ,其 序列 由 breakpoint 
之 间 的 外 显 子 的 碱 基 构 成 ;(2) 单一 外 显 子 , 即 
circRNA 位 于 已 知 基 因 的 某 一 个 外 显 子 内 ,其 序列 
由 breakpoint 之 间 的 所 有 碱 基 构 成 。(3 ) 基因 反 义 
链 , 即 circRNA 位 于 已 知 基因 的 反 义 链 ,其 序列 由 反 
义 链 上 breakpoint 之 间 的 所 有 碱 基 构 成 ;(4) 内 含 
子 , 即 circRNA 位 于 已 知 基因 的 某 一 个 内 含 子 内 ,其 











































































































的 柱状 图 。 
1.4 CircRNA 的 来 源 基 因 分 析 、 靶 向 预测 及 调控 
网 络 构 建 
通过 BLAST 将 circRNA 的 来 源 基 因 注 释 到 GO 
All KEGG 数据 库 ,从 而 对 cireRNA 的 功能 进行 预测 。 
为 研究 cireRNA 的 竞争 性 内 源 RNA ( competing 
endogenous RNA, ceRNA ) 作 用 ,利用 TargetFinder 软 
件 ( Allen et al., 2005) 预测 circRNA #8 [8] 25-4 AY 
miRNA ,得 到 circRNA 与 miRNA 之 间 的 靶 向 结合 3 
系 。 从 结果 中 提取 出 cireRNA 潜在 的 靶 向 结合 
miRNA 的 位 点 信息 并 按照 超 几 何 分 布 的 P 值 排序 ， 
选择 P 值 最 小 的 前 300 对 circRNA-miRNA 关系 , 利 
用 Cytoscape v.3.2.1 软件 (Odelius et al., 2011) #4 
建 调控 网 络 图 。 

进一步 利用 TargetFinder 软件 预测 circRNA 38 
向 miRNA 结合 的 mRNA ,根据 circRNA 与 miRNA, 
miRNA 与 mRNA 的 结合 关系 , 得 到 circRNA- 
miRNA-mRNA 的 靶 向 结合 关系 ,最 后 通过 Cytoscape 
v.3.2.1 软件 对 调控 网 络 进行 可 视 化 。 
1.5 CircRNA 的 RT-PCR 验证 

为 验证 预测 出 的 意 蜂 cireRNA ,随机 选取 3 个 
circRNA ( novel. cire; 000663, novel. cire, 002507 及 
novel cire, 004725) 进行 RT-PCR WIE, H Guo 等 
(2018) 的 方法 ,利用 DNAMAN 软件 根据 所 选 
circRNA 的 碱 基 序列 设计 相应 的 特异 性 背靠背 引物 
与 线性 扩 增 引物 ,委托 上 海 生 工 生 物 工程 有 限 公司 
进行 合成 ,引物 序列 见 表 1。 利用 RNA 抽 提 试剂 盒 
(TaKaRa 公司 ,中 国 ) 分 别提 取 Am7 和 Am10 样品 
的 总 RNA ,分 为 2 份 ,其 中 一 份 总 RNA 用 3 U/mg 
RNase R( 吉 赛 生 物 公司 ,中 国 ) 进行 消化 以 去 除 线 
性 RNA, 37% 处 理 15 min 后 以 随机 引物 进行 反 转 
录 , 得 到 circRNA 的 cDNA; 男 一 份 总 RNA 直接 作为 
模板 ,以 Oligo (dT)»w 引 物 进行 反 转 录 , 得 到 线性 
RNA 的 cDNA。 利 用 基因 组 DNA ( gDNA ) 抽 提 试 剂 

































































序列 由 breakpoint 之 间 的 所 有 碱 基 构 成 ;(5 ) 外 显 子 
和 内 含 子 , 即 cireRNA 的 breakpoint 之 间 有 已 知 基因 
的 一 个 或 多 个 完整 外 显 子 但 breakpoint 不 位 于 已 知 
外 显 子 的 起 始 或 终止 ,其 序列 由 上 breakpoint 起 始 
位 点 到 第 一 个 外 显 子 起 始 位 点 .breakpoint 之 间 的 外 
显 子 和 最 后 一 个 外 显 子 到 breakpoint 终止 位 点 的 碱 
基 构 成 ;(6) 基 因 间 区 , 即 cireRNA 位 于 已 知 基 因 的 
之 间 ,其 序列 由 breakpoint 之 间 的 所 有 碱 基 构 成 ; 对 
circRNA 的 来 源 基 因 在 西方 蜜蜂 染色 体 上 的 位 置 进 
行 统计 ,绘制 每 条 染色 体 的 cirecRNA 来 源 基因 数量 




































































f ( Axygen 公司 ,中 国 ) 分 别提 取 Am7 和 Am10 样品 
的 gDNA。 获 得 的 Am7 和 AmlO 的 cDNA 及 gDNA 
等 摩尔 混合 ,作为 PCR 的 模板 。PCR 体系 为 20 uL: 
模板 1 pL, Mixture 10 pL， 上 下 游 引 物 (10 pmol/ 
L) 各 1 pL, 无 菌 水 补 至 20 uL, PCR 反应 程序 设 
置 : 94%C 预 变性 5 min; 94C 变 性 50 s, 按 照 引 物 Tm 
值 对 应 的 退火 温度 30 s, 72% 延伸 1 min, 共 36 个 循 
环 ; 72°C 延伸 5 min, PCR 产物 经 1.5% 琼 脂 糖 凝 胶 
电泳 检测 。 
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表 1 RT-PCR 验证 的 引物 信息 reads 平均 为 136 463 071 条 ( 表 2); 各 样品 的 平均 
Table 1 Primer information in RT-PCR validation Q20 和 Q30 分 别 为 98. 48% 和 95.51% ( X 25 " E 
E Rd d 述 结 果 表 明 本 研究 中 的 数据 质量 良好 ,可 用 于 下 一 
ire name rimers rimer sequences 步 分 析 。 各 样品 的 isch ea 数 均 在 136 779 122 
Novel_cire 000663 DI-F | GGAAAGGCTGTGTTCAAGG p PR j B 
DI-R — CCTCTTCGTTGTTCGCTT 条 及 以 上 , 比 对 上 西方 蜜蜂 参考 基因 组 的 anchor 数 
CIF CCATTGGAACCGATTATTGG 均 在 17 502 228 条 及 以 上 ( 表 3)。 
CI-R CAGATTCAGGTGTTATCCCTTG 2.2 意大利 密 蜂 工蜂 中 肠 的 circRNA 的 预测 及 结 
Novel cire. 002507 D2-F ATTTCCTTGGGCATAGCC 构 特 征 
D2-R CTCGGTCAAACCATACACC 在 Am7 和 AmlO 中 共 预 测 到 10 833 个 
C2-F CGCTATCTTATGCCTCCCT 
a ae T r Xe ov ANS d 
CR — CTCAATGATTTGCTCTTGGC circRNA ,二 者 的 共有 circ RNA 数 为 5 600 个 ,特有 
. NEN : 
Novel cir 004725 D3-F CTGGGTTACCGTCGTATCCT circRNA 数 分 别 为 1741 和 3 492 个 (图 1) ;circRNA 
D3-R — TAGCAGGCTGAAGACCTCTC 的 长 度 最 小 仅 15 nt ,最 长 可 达 98 429 nt, 主要 分 布 
C3-F  CTACAATCCTTATCACCACGG 1E 15 ~ 1 000 nt (70. 99%) d F(A 2: B), 
C3-R  TAATGGCGGTCTCAGCAT 


DI-F/DI-R, D2-F/D2-R, D3-F/D3-R 为 背靠背 引物 ; CI-F/CI-R, 
C2-F/C2-R, C3-F/C3-R 为 线性 扩 增 引物 。D1-F/D1-R,，D2-F/D2-R 
and D3-F/D3-R are divergent primers; CI-F/CI-R, C2-F/C2-R and 
C3-F/C3-R are convergent primers. 


2 结果 
2.1 高 通 量 测 序数 据 的 质 控 与 评估 

本 研究 的 意大利 蜜蜂 工蜂 中 肠 样 品 的 测序 平均 
得 到 140 926 764 条 raw reads, 经 过 滤 得 到 clean 























CircRNA 的 类 型 中 ,已 注释 的 外 显 子 circRNA 最 多 ， 
达 5 261 个 (48. 5696) , 其 次 为 反 义 链 circRNA 
(1715 个 , 15. 83% ) .单一 外 显 子 circRNA (1 628 
个 , 15.03% ) .外 显 子 和 内 含 子 circRNA(1 141 个 ， 
10.53% ) N&T circRNA(780 个 , 7.2096 ) 以 及 基 
因 间 circRNA(308 个 , 2.8496 ) (图 2: A) 。 进 一 步 
分 析 发 现 ,分 布 于 意大利 蜜蜂 1 号 染色 体 上 的 
circRNA 数量 最 多 , 达 2 083 个 (19.79% ) ,其 次 为 8 
号 染色 体 (1 941 个 , 18.4496 ) 和 2 号 染色 体 (1 102 
个 , 10.17% )( 图 2.: C)。 



































表 2 意大利 蜜蜂 工蜂 中 肠 样品 RNA-seq 数据 的 信息 统计 
Table2 Summary of RNA-seq data in the midgut of Apis mellifera ligustica workers 
SA st pg trado p 
C ee Tta e) en) 
Am7-1 160 844 082 155 961 044 (96.96) 22 780 265 834 (98.52) 22 103 303 455 (95.59) 
Am7-2 129 878 194 126 476 532 (97.38) 18 477 147 770 (98.52) 17 930 815 636 (95.60) 
Am7-3 113 683 898 110 932 376 (97.58) 16 176 905 690 (98.43) 15 679 232 245 (95.40) 
Aml0-1 160 537 248 155 007 614 (96.56) 22 626 100 922 (98.51) 21 960 409 097 (95.61) 
Am10-2 149 230 808 144 360 740 (96.74) 21 060 425 866 (98. 46) 20 421 732 789 (95.47) 
Am10-3 131 386 354 126 040 122 (95.93) 18 379 072 115 (98.44) 17 811 645 341 (95.40) 


Q20: 99% fd 5E TE WHR 99% base correct rate; Q30; 99.9% li XE TE WAR 99.9% base correct rate. 


X3 未 比 对 上 核糖 体 数据 库 的 clean reads 比 对 西方 蜜蜂 参考 基因 组 的 信息 统计 


Table3 Summary of the mapping information of unmapped clean reads in ribosome database in 





Apis mellifera reference genome 





样品 rRNA 去 除 率 (% ) 短 序列 读 自 比 对 上 的 短 序列 比 对 率 (% ) 
Samples rRNA removal rate Anchor reads Mapped anchors Mapping rate 
Am7-1 0.08 236 400 908 19 006 601 8.04 
Am7-2 0.09 168 420 472 19 873 468 11.80 
Am7-3 0.07 136 779 122 17 502 228 12.80 
Am10-1 0.14 195 020 440 22 441 139 11.51 
Aml10-2 0.11 184 091 234 21 109 367 11.47 
Aml10-3 0.12 168 790 774 17 765 330 10.53 
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3492 (32.23%) 


1741 (16.08%) 


5600 (51.69%) 





Am7 Am10 























A1 意大利 蜜蜂 7(Am7) 和 10 (Aml0) 日 龄 工蜂 
中 肠 circRNA 的 韦 恩 分 析 
Fig. 1 Venn diagram of circRNAs in the midguts of 


7-old-day ( Am7) and 10-old-day (Am10) adults 


of Apis mellifera ligustica workers 
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2 意大利 蜜蜂 工蜂 中 上 肠 circRNA 的 结构 特征 


Number 


Jt 
AL 


BL 








2.3 意大利 蜜蜂 工蜂 中 肠 circRNA 的 来 源 基 因 的 
GO 分 类 


GO 分 类 结果 显示 , 意大利 蜜蜂 工蜂 中 上 肠 
circRNA 的 来 源 基 因 可 注释 到 45 个 GO term ,涉及 
生物 学 进程 .分子 功 能 和 细胞 组 分 三 大 功能 分 类 
(图 3)。 生 物 学 进程 中 ,注释 基因 数 前 5 位 的 GO 
term 分 别 是 细胞 进程 (679 个 ) 代谢 进程 (633 个 )、 
单 组 织 进程 (547 个 )、 定 位 (195 个 ) 及 生物 调控 
(188 个 ) ;分 子 功能 中 ,注释 基因 数 前 5 位 的 分 别 是 
结合 (734 个 ) .催化 活性 (563 个 ) .转运 器 活性 (98 
个 ) 分 子 传感器 活性 (80 个 ) 及 核酸 结合 转录 因子 
活性 (67 个 ) ;细胞 组 分 中 ,注释 基因 数 前 5 位 的 分 
别 是 细胞 (334 A) ` 细 胞 组 分 (334 个 )、 细 胞 膜 (308 
个 ) 细胞 膜 组 分 (297 个 ) 及 细胞 器 (234 个 )( 图 3)。 








E 已 注释 外 显 子 circRNA Annotated _exonic circRNA 
m 单一 外 显 子 circRNA Single _exonic cireRNA 

Il 反 义 链 circRNA Antisense circRNA 

E 内 含 子 circRNA Intronic circRNA 

m 外 显 子 和 内 含 子 circRNA Exonic _intronic cireRNA 
m 基因 间 circRNA Intergenic circRNA 
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1 3 
西方 蜜蜂 染色 体 Chromosomes of Apis mellifera 





Fig. 2 Structural characteristics of circRNAs in the midgut of Apis mellifera ligustica workers 


A; CircRNA 的 类 型 Type of circRNAs; B: CircRNA 的 长 度 分 布 Length distribution of circRNAs; C; CircRNA 在 西方 蜜蜂 1 10 号 染色 体 上 的 分 


布 Distribution of circRNAs on chromosomes 1 - 10 of A. mellifera. 


上 述 结果 说 明 circRNA 广泛 参与 意 蜂 工蜂 中 肠 的 新 
陈 代谢 和 细胞 生命 活动 的 诸多 方面 。 
2.4 意大利 蜜蜂 工蜂 中 肠 circRNA 的 来 源 基 因 的 
KEGG 数据 库 注释 

进一步 对 意大利 蜜蜂 工蜂 中 肠 的 circRNA 的 来 
源 基因 进行 KEGG 数据 库 注释 ,结果 显示 涉及 121 
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条 代谢 通路 ,其 中 注释 基因 数 前 10 位 分 别 是 内 否 作 
用 (39 个 ) 内 质 网 蛋白 加 工 (37 个 )、 核 糖 体 (33 
个 ) Woa 信号 通路 (31 个 )、 味 叭 代谢 (30 个 ) 、 胰 岛 
素 耐 受 性 (29 个 ) 、 碳 代谢 (29 个 ) Hippo 信号 通路 
(28 个 ) ,FoxO 信号 通路 (23 个 ) 以 及 神经 活性 配 体 
受 体 交互 (21 个 )( 表 4) ,再 次 说 明 circRNA 广泛 参 
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图 3 意大利 蜜蜂 工蜂 中 肠 circRNA 的 来 源 基因 的 GO 数据 库 注 释 


Fig. 3 GO database annotations of the source genes of circRNAs in the midgut of Apis mellifera ligustica workers 
与 意 蜂 中 肠 的 生长 发育 和 代谢 等 生物 学 过 程 。 靶 向 结合 127 个 miRNA, ,其 中 绝 大 多 数 (70. 65% ) 
2.5 意大利 蜜蜂 工蜂 中 肠 circRNA-miRNA 调控 仅 能 结合 1 个 miRNA, tF, novel_ circ _003007, 
网 络 的 构建 及 分 析 novel_circ_003010, ，novel_circ_004674 , novel_cire_ 
软件 预测 结果 显示 ,共有 1 707 个 circRNA 可 ”013072 及 novel_cire_013073 结合 的 miRNA 数 最 多 ， 
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表 4 意大利 蜜蜂 工蜂 中 肠 circRNA 的 来 源 基因 注释 数 前 20 位 的 KEGG 代谢 通路 


Table 4 


Top 20 KEGG pathways annotated by circRNAs' source genes in the midgut of 


Apis mellifera ligustica workers 
































通路 基因 数 @ 值 
Pathway Gene number Q value 
内 知 作 用 Endocytosis 39 0.164158 
内 质 网 蛋白 加 工 Protein processing in endoplasmic reticulum 37 0. 444072 
核糖 体 Ribosome 33 0.700042 
Wat 信号 通路 Wnt signaling pathway 31 0.001494 
IE] Purine metabolism 30 0. 968456 
胰岛 素 耐 受 性 Insulin resistance 29 0. 000003 
碳 代谢 Carbon metabolism 29 0. 700042 
Hippo 信和 号 通路 Hippo signaling pathway 28 0. 000746 
FoxO 信号 通路 FoxO signaling pathway 23 0.052567 
神经 活性 配 体 受 体 交 互 Neuroactive ligand-receptor interaction 21 0. 107067 
吞噬 体 Phagosome 21 0.348897 
磷脂 酰 肌 醇 信号 系统 Phosphatidylinositol signaling system 20 0. 094723 
氨基 酸 生 物 合成 Biosynthesis of amino acids 19 0.533168 
溶 酶 体 Lysosome 19 0.700042 
泛 素 介 导 蛋白 水 解 Ubiquitin-mediated proteolysis 19 0. 999999 
甘油 磷脂 代谢 Glycerophospholipid metabolism 17 0.700042 
TGF-beta 信号 通路 TGF-beta signaling pathway 16 0.060735 
寿命 调节 通路 Longevity regulating pathway 16 0. 345036 
糖 酵 解 / 糖 质 新 生 Glycolysis/Gluconeogenesis 16 0.402383 
mRNA 监督 通路 mRNA surveillance pathway 16 0.935653 


均 为 8 个 ;此 外 ,novel_circ_013374, novel_circ_ 
013172, novel. circ, 012110, novel, circ, 011938, 
novel circ, 011834 和 novel. circ, 011770 可 分 别 结合 
7,6,5,5,4 和 2 个 miRNA( 图 4)。 上 述 结果 说 明 
意大利 蜜蜂 工蜂 中 肠 circRNA 可 通过 吸附 结合 
miRNA 对 下 游 的 基因 表达 产生 调控 作用 。 
2.6 意大利 蜜蜂 工蜂 中 肠 circRNA-miRNA- 
mRNA 调控 网 络 构 建 及 分 析 
进一步 构建 意大利 蜜蜂 工蜂 中 肠 circRNA- 
miRNA-mRNA 调控 网 络 , 结 果 显 示 141 个 circRNA 
与 110 个 miRNA 和 283 个 mRNA 之 间 形 成 两 个 复 
杂 的 调控 网 络 ,多 数 miRNA 处 于 调控 网 络 的 内 部 ， 
而 circRNA 与 mRNA 多 处 于 调控 网 络 的 外 周 ， 
miRNA 可 同时 结合 多 个 circRNA 和 mRNA (KI 5), 
再 次 说 明 意大利 蜜蜂 工蜂 中 肠 cire RNA 通过 竞争 性 
地 结合 miRNA 影响 miRNA 与 mRNA 的 结合 ,进而 
对 下 游 的 基因 表达 产生 调控 。 

对 上 述 调控 网 络 相 关 的 mRNA 进行 GO 分 类 ， 
结果 显示 它们 共 涉 及 28 个 GO term, 包 括 细胞 进程 
(86 个 ) 和 代谢 进程 (68 个 ) 等 10 个 生物 学 进程 相 
XK term ,结合 (95 个 ) 和 催化 活性 (54 个 ) 等 9 个 分 






















































































子 功能 相关 term ,细胞 (33 个 ) 和 细胞 组 件 (33 个 ) 
等 9 个 细胞 组 分 相关 term。 进 一 步 的 KECC 代谢 通 
路 富 集 分 析 结果 显示 ,上 述 mRNA 涉及 150 条 代谢 
通路 ,包括 鞘 糖 脂 (8 个 ) 和 粘 多 糖 (4 个 ) 等 27 条 新 
陈 代谢 相关 通路 ,多 巴 胺 突 触 (4 个 ) 和 胰岛 素 信 和 号 
通路 (3 个) 等 47 条 有 机 体系 统 相关 通路 ,mRNA 监 
视 通 路 (3 个 ) 和 剪接 体 (3 个 ) 等 7 条 遗传 信息 加 工 
相关 通路 ,cAMP 信号 通路 (8 个 ) 和 Hippo 信号 通路 
(5 个 ) 等 24 条 环境 信息 加 工 相 关 通 路 ,以 及 溶 酶 体 
(5 个 ) 和 细胞 凋 亡 (1 个 ) 等 9 条 细胞 进程 相关 通 
路 。 上 述 结 果 说 明 circRNA 可 通过 ceRNA 网 络 在 
意 蜂 工蜂 中 肠 的 物质 代谢 、 能 量 代谢 .生长 发 育 、 信 
号 通路 和 免疫 防御 的 诸多 方面 发 挥 作用 。 
2.7 意大利 蜜蜂 工蜂 中 上 肠 circRNA 的 RT-PCR 验证 
随机 挑选 的 3 个 circRNA ( novel. cire. 000663 , 
novel. circ, 002507 及 novel_cire_004725 ) ,设计 对 应 
的 背靠背 引物 和 线性 扩 增 引物 (图 6: A) 进 行 RT- 
PCR 验证 。 电 泳 结果 显示 ,以 cDNA 作为 模板 ,靠背 
引物 和 线性 扩 增 引物 均 能 扩 增 出 目的 片段 ; 以 
gDNA 作为 模板 ,线性 扩 增 引物 能 够 扩 增 出 目的 片 
段 ,但 背靠背 引物 不 能 扩 增 出 目的 片段 (图 6: B- 





















































































































































































































































12 期 


Eg 


EUST AE. 意大利 蜜蜂 工蜂 中 肠 的 环 状 RNA 及 其 调控 网 络 分 析 1371 





ant 
We 











" 
ae a 
uy L2 
Shia 011834 
eem rag; novel $ire/003010 


M, lax 0300 ~ 
XXX | 
"ea a E 

; ur 


c IIIIIIIIl 


novel 02580 


图 4 意大利 蜜蜂 工蜂 中 肠 circRNA-miRNA 调控 网 络 
Fig. 4 Regulatory networks between circRNAs and miRNAs in the midgut of Apis mellifera ligustica workers 
黄色 椭圆 代表 circRNA ;绿色 三 角形 代表 miRNA, Yellow ellipses indicate circRNAs, and green triangles indicate miRNAs. 





图 5 意大利 蜜蜂 工蜂 中 肠 circRNA-miRNA-mRNA 调控 网 络 
Fig. 5 Regulatory networks between circRNAs, miRNAs and mRNAs in the midgut of Apis mellifera ligustica workers 
黄色 椭圆 代表 circRNA ;绿色 三 角形 代表 miRNA; 20 (68 2:JE [XX mRNA, Yellow ellipses indicate circRNAs, green triangles indicate miRNAs, and 
red diamonds indicate mRNAs. 
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D) ,一 定 程 度 证 明 本 研究 预测 出 的 cireRNA 真实 
存在 。 
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图 6 意大利 蜜蜂 工蜂 中 肠 cireRNA AY RT-PCR 验证 
Fig. 6 RT-PCR validation of circRNAs in the midgut 








of Apis mellifera ligustica workers 

A: 背靠背 引物 和 线性 扩 增 引物 的 扩 增 原理 示意 图 Diagram of 
amplification with divergent and convergent primers; B - D; 分 别 是 
Novel_circ_000663 ，Novel_circ_002507 和 Novel_circ_004725 的 背 靠 
背 引物 和 线性 扩 增 引物 扩 增 产物 的 琼脂 糖 凝 胶 电 泳 Agarose gel 
electrophoresis for amplification products with divergent and convergent 
primers of Novel_cire_000663 , Novel_cire_002507, and Novel_cire_ 
004725 , respectively. 





3 讨论 


目前 ,circRNA 的 研究 主要 集中 在 人 类 (Jeck et 
al., 2013; Salzman et al., 2013; Memczak et al., 
2013) .哺乳 动物 (Gruner et al., 2016; Liang et al., 
2017) 和 重要 疾病 (Spitali and Aartsma-Rus, 2012; 
Taulli et al., 2013; Hsiao et al., 2017; Zhu et al., 
2017) 等 领域 。 随 着 高 通 量 测序 技术 和 生物 信息 学 
分 析 方 法 的 高 速 发 展 和 应 用 ,人 们 已 在 如 水 稻 (Lu 
et al., 2015) R tẹ ( Westholm et al., 2014) , 3% 4 
( Gan et al., 2017) 、 古 细菌 (Danan et al., 2012) 和 球 
FEW (Guo et al., 2018) 等 越 来 越 多 物种 中 预测 和 鉴 








XE Hj circRNA, Westholm 等 (2014) Xf Œ tE AJ 
circRNA 进行 全 基因 组 分 析 , 检测 到 了 2 500 个 
circRNA ,大 部 分 circRNA 由 外 显 子 环 化 而 来 。GCan 
等 (2017) 通 过 对 家 看 丝 腺 进行 高 通 量 测序 , 共 检 测 
到 3 155 个 circRNA ,其 中 50.26% 为 外 显 子 和 内 会 
子 环 化 circRNA, Chen 等 (2018 ) 通过 高 通 量 测序 
技术 和 生物 信息 学 方法 从 意大利 蜜蜂 蜂王 的 卵 梨 组 
织 中 预测 出 12 211 个 cireRNA ,超过 80% 由 外 显 子 
环 化 ,这 是 迄今 唯一 有 关 蜜 蜂 circRNA 的 研究 报道 。 
物种 的 基因 组 中 仅 有 部 分 基因 能 够 转录 形成 
circRNA( Jeck and Sharpless, 2014; Ebbesen et al., 
2016) 。 对 于 人 类 ,大 约 有 5.8% ~ 23% 的 基因 在 不 
同 的 组 织 和 细胞 类 型 中 转录 形成 cricRNA( Conn et 
al., 2015; Kelly et al., 2015 )。 在 猪 胚 发 育 过 程 中 
的 不 同时 期 的 大 脑 样品 中 转录 形成 circRNA 的 基因 
比例 在 5.8% ~15.8 之 间 (Morten et al., 2015)。 本 
研究 利用 circRNA-seq 技术 对 意大利 蜜蜂 7 和 10 日 
龄 工蜂 中 肠 进行 测序 , 共 预 测 出 10 833 个 circRNA, 
其 中 由 已 注释 的 外 显 子 环 化 而 来 的 circRNA 数量 最 
多 ,达到 5 261 个 ,与 果 晶 (Westholm et al., 2014) 及 
意大利 蜜蜂 蜂王 circRNA( Chen et al., 2018) 的 研究 
结果 一 致 , 但 与 家 看 circRNA( Gan et al., 2017) 的 研 
究 结果 存在 差异 。 各 样品 的 anchor reads 比 对 上 西 
方 蜜蜂 参考 基因 组 ( Honeybee Genome Sequencing 
Consortium, 2006) 的 比例 介 于 8. 04% ~ 12. 8096 
( 表 3) ,与 其 他 物种 的 研究 结果 类 似 。 此 外 ,大 部 分 
anchor reads 未 比 对 上 参考 基因 组 ,表明 Bowtie2 软 
件 (Langdon, 2015) 比 对 参考 基因 组 的 筛选 标准 及 
西方 密 蜂 基因 组 的 组 装 仍 有 待 优 化 。 本 研究 预测 出 
的 circRNA 可 为 意大利 蜜蜂 的 ncRNA 信息 提供 重 
要 补充 ,也 为 其 他 蜂 种 的 cireRNA 研究 提供 了 有 益 
的 参考 信息 。 

CircRNA 的 表达 具有 发 育 时 期 和 组 织 特异 性 
(Li et al., 2017) 。 本 研究 中 ,共有 5 600 个 circRNA 
在 7 和 10 日 龄 成 年 工蜂 中 肠 中 共同 表达 , 特 
达 的 circRNA 数 分 别 为 1 741 和 3 492 个 ,推测 这 
共同 circRNA 在 意大利 蜜蜂 工蜂 中 肠 发 育 ies 
挥 基础 性 的 调节 功能 ,而 特有 circRNA 在 中 肠 发 育 
的 不 同 阶段 发 挥 特殊 的 调控 作用 。 然 而 ,本 研究 预 
测 出 的 cireRNA 是 否 在 其 他 日 龄 的 意大利 蜜蜂 工蜂 
中 肠 表 达 , 以 及 是 否 在 意大利 蜜蜂 工蜂 的 其 他 组 织 
和 上 需 官 表达 仍 需 要 进一步 研究 。 

本 研究 中 ,G0 分 析 绪 果 显示 ,分 别 有 633, 12 
和 3 个 cireRNA 的 来 源 基因 涉及 意大利 蜜蜂 工蜂 中 
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肠 的 代谢 进程 .发 育 进 程 及 生长 (图 3) ,表明 相应 的 
circRNA 参与 调控 中 上 肠 的 生长 ,发育 和 新 陈 代谢 ;分 


工蜂 中 上 肠 circRNA-miRNA -mRNA 调控 网 络 ,分 析 结 
果 显 示 三 者 之 间 存 在 错综复杂 的 调控 关系 ,circRNA 


























BNA 151 和 3 个 cireRNA 来源 基因 涉及 应 激 反应 和 
免疫 系统 进程 ,表明 相应 的 circRNA 在 中 肠 的 免疫 
防御 中 发 挥 调 控 功 能 。KEGG 代谢 通路 富 集 分 析 结 
果 显 示 , citcRNA 的 来 源 基 因 参 与 意大利 蜜蜂 工蜂 
中 肠 的 121 条 代谢 通路 ,包括 氨基 酸 生 物 合成 (19 
个 基因 ) 和 mRNA 监视 通路 (16 个 基因 ) 等 物质 代 
谢 通 路 , 碳 代 谢 (29 个 基因 ) 和 寿命 调节 通路 (16 个 
基因 ) 等 能 量 代谢 通路 , 内 否 作 用 (39 个 基因 ) 和 
FoxO 信和 号 通路 (23 个 基因 ) 等 免疫 通路 ( 表 4)。 
Hippo 信号 通路 在 很 多 组 织 中 通过 促进 细胞 增殖 和 
抑制 细胞 凋 亡 来 调 探 需 官 大 小 ,在 发 育 中 发 挥 重 要 
作用 (Zhao et al., 2011) ; Wot 信号 通路 在 器 官 形成 
和 细胞 更 新 中 扮演 着 重要 角色 (Jansson et al., 
2015 ) 。 本 研究 发 现 ,意大利 蜜蜂 工蜂 中 肠 circRNA 
的 来 源 基因 中 ,分 别 有 28 和 31 个 富 集 在 Hippo 和 
Wot 信号 通路 ,表明 相应 的 cireRNA 通过 调控 此 
2 条 信和 号 通路 对 意 蜂 工 蜂 中 肠 的 生长 和 发 育 产 生 
影响 。 

Salmena 等 (2011 ) 提出 了 ”竞争 性 内 源 RNA” 
假说 , 即 任 何 包含 miRNA 反应 元 件 (miRNA 
response element) 的 RNA ,包括 mRNA、 假 基因 、 
IncRNA 和 circRNA ,可 以 竞争 性 地 结合 miRNA。 此 
后 , 越 来 越 多 的 研究 结果 都 证 实 了 该 假说 (Deng et 
al., 2016; Liu et al., 2016; Pan et al., 2017; Xu et 
al., 2017) , Pan 5$ (2017 ) 通过 对 冠 心病 人 和 正常 
人 中 DEcircRNA. 进行 分 析 , 发现 9 个 差异 表达 的 
circRNA n] DL 3B ix: 8B [8] Zi g hsa-miR-130a-3p ,从 而 
促进 TRPM3 ( transient receptor potential cation 
channel subfamily M member 3) 的 表达 和 参与 修复 血 
管 的 调控 。Liu 等 (2016) 通 过 研究 与 骨 关 节 炎 相关 
的 circRNA , 发 现 circRNA-CER 可 以 充当 miR-136 
的 “海绵 ” ,通过 吸附 后 者 对 基质 金属 蛋白 酶 13 编 
人 码 基 | 因 ( matrix. metallo proteinase 13 gene, MMP13) 
的 表达 进行 调控 而 参与 软骨 细胞 外 基质 的 降解 。 本 
研究 中 ,多 达 1 707 个 意大利 蜜蜂 工蜂 中 肠 circRNA 
up Sg] ii 127 个 miRNA ,其 中 部 分 circRNA n] $8 
向 多 个 miRNA ,如 novel_circ_003007 及 novel_cire_ 
003010 可 结合 8 个 miRNA; 也 有 多 个 circRNA 共同 
丢 向 同一 个 miRNA, 如 novel. cire. 010951, novel . 
circ, 002283 和 novel. cire, 002616 共同 靶 向 mir- 
8485 -y ;novel cire, 002579 和 novel, cire, 002580 共 
同 革 向 ame-miR-6066。 进 一 步 构建 的 意大利 蜜蜂 













































































































































































































































































































































































可 与 mRNA 竞争 性 地 结合 miRNA。 通 过 GO 分 类 
F KEGG 代谢 通路 富 集 分 析 进 一 步 预测 意大利 密 
蜂 工 蜂 中 肠 circRNA 的 功能 ,结果 显示 调控 网 络 中 
的 mRNA 涉及 细胞 进程 (86 个 ) 代谢 进程 (68 个 )、 
生物 学 调控 (25 个 ) 应 激 反应 (18 个 )、 信 号 通路 
(16 个 ) 和 发 育 进程 (2 个 ) 等 28 个 GO term; 此 外 
KEGG 代谢 通路 富 集 分 析 结 果 显 示 , 还 涉及 鞘 糖 脂 
(8 个 ) 粘 多 糖 (4 个 ) ,多 巴 胺 突 触 (4 个 ) .剪接 体 
(3 个 ) 溶 酶 体 (5 个 ) ZUR T7 CY 个 ) 内 吞 作用 
(1 47) , Wnt 信号 通路 (4 个 ) Hippo 信号 通路 (2 
个 ) MAPK 信号 通路 (2 个 ) 和 NF-kappa B 信和 号 通 
路 (1 个 ) 等 150 条 代谢 通路 ;上 述 结 果 暗 示 circRNA 
不 仅 能 通过 作用 于 来 源 基因 也 能 够 通过 作为 
ceRNA 影响 意大利 蜜蜂 工蜂 中 肠 的 生长 发 育 和 免 
疫 防 御 。 

CircRNA 除 可 直接 调控 来 源 基 因 (Zhang et al., 
2013; Li et al., 2015) 的 表达 外 ,还 能 够 调控 基因 的 
n] Æ BY E ( Ashwal-Fluss et al., 2014) 和 转录 过 程 
(Zhang et al., 2014) , 以 及 吸附 核糖 体 结 合 蛋 白 
(ribosome binding proteins, RBPs) ( Ashwal-Fluss et 
al., 2014) M/E RBPs 复合 体 的 支架 ( Qu et al., 
2015) 发 挥 作 用 。 最 新 研究 发 现 ,一 些 circRNA 具有 
编码 蛋白 的 能 力 (Legnini et al., 2017)。 推 测 本 人 研 
究 预 测 出 的 circRNA 也 能 通过 多 种 途径 参与 对 意 大 
利 蜜蜂 工蜂 中 上 肠 的 生长 和 发 育 过 程 的 调控 ,其 具体 
功能 有 待 于 进一步 深入 研究 。 

综 上 所 述 ,本 研究 结合 高 通 量 测序 技术 和 生物 
言 息 学 分 析 方 法 对 意大利 蜜蜂 工蜂 中 肠 circRNA git 
行 预测 和 分 析 ,研究 结果 提供 了 中 肠 cireRNA 的 数 
量 种 类 结构 特征 和 作用 的 信息 ,为 蜜蜂 的 ncRNA 
信息 提供 了 有 益 补 充 ,也 为 深入 研究 circRNA 在 中 
肠 发 育 及 胁迫 响应 过 程 的 作用 建立 基础 。 下 一 步 将 
对 不 同日 龄 的 意大利 蜜蜂 工蜂 中 肠 circRNA 进行 差 
异 表达 分 析 和 调控 网 络 分 析 , 从 而 筛选 出 与 中 肠 生 
长 发 育 密切 相关 的 DEcircRNA 。 
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